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ABSTRACT 

We present results of a numerical simulation of the 3-state Potts model 
with random bond, in two dimension. In particular, we measure the critical 
exponent associated to the magnetization and the specihc heat. We also 
compare these exponents with recent analytical computations. 
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These last years, many studies have been devoted to the problem of the effect of 
randomness on 2 dimensional statistical models. The effect of randomness is supposed to 
be directly related to the critical exponent of the specihc heat, a, according to the well 
known Harris criterion |I|]. If a is positive then the disorder will be relevant, i.e. under 
the effect of the disorder, the model will reach a new critical behavior at a new critical 
point. Otherwise, if a is negative, disorder is irrelevant, the critical behavior will not 
change. The Harris criterion does not give any information on the behavior of models 
with a = 0 (marginal case), and the most studied case, the 2D Ising model, falls in this 
category. Due to its apparent simplicity, this model is the most studied example of a 
disordered model with many analytical results i B i B B 0 as well as numerical ones 

H, a m o. 



In the present work, we will focus on the 3-state Potts model with disorder. The 
pure 3-state Potts model having a positive a (a = 1/9), it is expected that a new 
critical behavior will be obtained. The main purpose of the present work is to study 
numerically this new behavior. A certain number of analytical studies of this model 
predict new critical exponents [^, ^ |^. In these papers, the 3-state Potts model is 

considered like a generalization of the Ising model by shifting a regularization parameter 
(namely the central charge of the corresponding conformal field theory describing the 
behavior of the pure model at the critical point.) Computations are performed by using a 
perturbed conformal held theory with a near marginal operator associated to the disorder. 
Predictions for critical exponents of the 3-state Potts model in presence of disorder are 
obtained by employing perturbative expansion with the machinery of the renormalisation 
group. In particular, the exponent associated to the spin-spin correlation function is given 
by 



with 



< a(0)a(/2) >~ 



2A' = 2A. 
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+ 0 ( 6 ^), 
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where e is the perturbation parameter. It measures the deviation from the Ising model, 
see O, Q for details.) The critical exponent associated to the energy-energy correlation 
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function is given by 



(3) 



< e:(0)£(i?) >~ R 



-2A, 



with 



2A^ = 2A^ - 3e + + O(e^). 



Here, e takes the value 



e = 



15 



(4) 



(5) 



for the 3-state Potts model. Plugging this value into the previous expressions, we obtain 



A' = — + 0, 00132 + O(e^) = 0, 13465 + O(e^). 
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It will be more convenient to express A'^ in function of A^ = Then 



( 6 ) 



A: = A.(1 + <5) 



6 = 0.0099 ~ 0.01 . 



For Ag, we have 



A' = l,02 + O(e^) 

which can be compared to the value for the pure 3-state Potts model 

A, = 0.8 . 



(7) 



( 8 ) 



(9) 



It is clear at this stage that the deviation of A^- due to the disorder will be very difficult 
to measure, while the one of A^ should be observed easily. 

A^ and A^ do also appear in the magnetization and in the specihc heat. Using some 
standard finite size scaling arguments [^, we have at the critical point 



M{L) ~ L~^- 
C{L) ~ 



( 10 ) 

( 11 ) 



These results were in fact obtained by using perturbative computations around the 
replica symmetry solntion. It can also be argued that it is necessary to break the replica 
symmetry, in a fashion similar to the Paris! solution for the mean field spin glass [|^ . 



If snch a solntion is employed, then the modihcation appears only at the third order in 



3 



e for the energy-energy correlation function, while there is no modihcation for the spin- 
spin correlation function (at least up to fourth order in e.) The new exponent for the 
energy-energy correlation function is given by |T5 



A: = A,-5£ + 0(£») = 1 + 0(£=). (12) 

Thus we have a certain number of exponents that we can compare to the numerical 
results that will be presented below. 

The Hamiltonian of the simulated model is given by 

« = - E (13) 

{iJ} 

where the coupling constant between nearest neighbor spins takes the value 



Jij pJo T (1 p)Ji- 



(14) 



Without any lost of generality, we can consider the case where p = \- Then the model is 
self-dual and thus the critical temperature is exactly known. It is given by the solution 
of the equation 



1 — e 

1 -I- (g — l)e~Wo 



= 



(15) 



In this Letter, we will only consider the case with a strong disorder, he. Jo = 1, Ji = 
The reason for such a choice is that there exists a cross-over length /j, depending on 
the strength of the disorder. It is only for distance larger than Ij that we expect to 
measure the critical behavior of the model with disorder. Extrapolating the results of 
the 2d random bond Ising model, we expect that with the present disorder, the cross-over 
length will take the value /j ~ 2 — 5. We will come back later to the manifestation of 
this cross-over length. Thus by running on large lattice, up to 1000^ for results presented 
here, we are quite sure to measure the critical behavior of the model with disorder. 



Monte Carlo data were obtained by using the well known Wolff cluster algorithm |]^ 



as well as the Swendsen- Wang algorithm . Details on the algorithms and the simula- 
tion parameters will be presented in a subsequent paper. Measurements were performed 
on a square lattice with helical boundary conditions. Due to the very strong disorder 
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that we consider, we needed to have huge statistics over the number of configurations 
of disorder. Simulations were performed for lattice with size ranging from L = 2 to 
L = 1000. For the magnetization, the number of conhgurations of disorder were 100000 
for L = 2 — 50, then 40000 for L = 100 to 4000 for L = 1000. For each of these con- 
hguration of disorder, measurements were taken over ti = 1000 conhgurations. For the 
specihc heat, the number of conhgurations need to be larger, namely t 2 = 20000. 10000 
conhgurations of disorder were simulated for L = 2 — 50, then from 1000 for L = 100 
to 100 for L = 1000. Finally, statistical errors were computed by taking the mean value 
over the conhgurations of disorder. The values of t\ and ^2 were determined in order 
that the thermal huctuations are smaller than the one coming from the distribution of 
disorder. 

The hrst result that we present is the Log-Log plot of the magnetization versus the 
lattice size L at the critical point (see Fig. 1). This plot exhibits a perfect scaling behavior 
with a very small deviation between L = 2 and L = 5. It is the manifestation of the 
cross-over length and justihes our previous assumption on its value. This is even more 
obvious on the second hgure. There we plot vs. ln(L). We see a jump in the 

cross-over region for 2 < L < 5. Outside this region we observe a rather good plateau. 
We now compute A'^ inside the plateau region with the parametrization of Eq. : 



At this level of precision, it is too difficult to draw any conclusion on the validity of the 
computed value vs. the measured one. The only obvious fact is that A'^ is very close to 
Ag-. With the statistics that we present here, it is not possible to differentiate these two 
exponents. 

For the specihc heat, the situation is different and we see clearly the inhuence of 
disorder. For the pure 3-state Potts model 



A'^ = 0.1337 ± 0.0007 = A^(l + 5) 



h = 0.003 ±0.005 . (16) 



This is to be compared with the analytical value 



6 = 0.01 . 



( 17 ) 




(18) 
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while for the model with disorder we obtain a negative exponent (see Fig. 3). The only 
possible parametrization of the cnrve is the one with a cnsp, i.e. 

C{L) =A + BL~\ (19) 

With snch a parametrization, we obtain 

a: = 0.13 ±0.04 (20) 

Using the relation x = 2A(, — 2 gives 

= 1.065 ±0.02 . (21) 

This resnlt shonld be compared with the analytical valnes 

A( = 1.02±O(e^) and A" = l±0(e^). (22) 



In conclnsion, we have measnred the critical exponent of the magnetization of the 3 
state Potts model in the presence of disorder. We have shown that its valne does not 
differ signihcantly from the pnre case valne. Comparisons with analytical compntations 
of this exponent [^, are beyond the level of precision of the present simnlation. For 
the critical exponent associated to the specihc heat (or the energy-energy correlation 



fnnction), the nnmerical valne obtained in Eq. is very close to the analytical one ||^ 



0 . 0 . 13 - The resnlt given by the nnmerical simnlation seems to agree better with the 
analytical valne obtained with a replica symmetry ansatz (A( = 1.02 ± 0{e^)) than the 
one given by the replica symmetry breaking ansatz A" = 1 ± O(e^). 
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Note 



As this work was being completed I became aware of a related work by J.-K. Kim 
which deals with a similar problem. In particnlar, Kim also measnred the exponent of 
the magnetisation (the parameter rj of the magnetic snsceptibility in fact) for the case 
with weak disorder and hnds it nnchanged from the case withont disorder. 
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Figure Captation 



Fig. 1 \nMag{L) vs. InL for 3-state Potts model with disorder 1/10 
Fig. 2 Mag{L) x vs. InL for 3-state Potts model with disorder 1/10 
Fig. 3 C{L) vs. L for 3-state Potts model with disorder 1/10, with a best fit. 
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